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Abstract 


The Weizsacker- Williams (WW) method of virtual quanta is used to make 
approximate cross section calculations for peripheral relativistic heavy-ion col- 
lisions. We calculated the Coulomb fission cross sections for projectile ions of 
208 Pb and 209 Bi with energies of 158 A GeV interacting with a 208 Pb target. 

We also calculated the electromagnetic absorption cross section for a 208 Pb ion 
interacting as described. For comparison we use both the full WW method and 
a standard approximate WW method. The approximate WW method resulted 
in larger cross sections compared with the more accurate full WW method. 

1. Introduction 

An important process neglected in the transport codes FLUKA (refs. 1 and 2) and 
MCNPX (ref. 3) is nuclear electromagnetic dissociation (refs. 4-8). When two nuclei 
collide, with an impact parameter less than or equal to the sum of the radii, they break 
up due to the strong forces. However, if the impact parameter is greater than the sum of 
the nuclear radii then breakup can occur via the electromagnetic (EM) interaction. This is 
especially important for few-nucleon removal (including neutrons) and for medium nuclei 
such as Aluminum or Iron. The few-nucleon EM removal cross sections can be larger 
than strong interactions cross sections (ref. 5). The EM interaction also leads to large 
fission cross sections (refs. 6 and 7) and double EM processes lead to a copious amount of 
electron-positron production (ref. 4) with cross sections in the kilobarn region. FLUKA 
and MCNPX include none of these processes. 

The transport code HZETRN uses the nuclear fragmentation code NUCFRG (ref. 
9), which contains a description of high energy nucleus-nucleus collisions in terms of a 
modified abrasion- ablation model. An advantage of NUCFRG is that it contains code for 
nuclear break up due to both the strong and electromagnetic interactions (ref. 8). 

None of the above three transport codes include electromagnetic dissociation cross 
sections leading to few-nucleon removal or fission. The present paper continues the study 
of relativistic Coulomb fission (refs. 6 and 7) with a view to including Coulomb fission 
cross sections in future versions of transport codes. 

There are two basic mechanisms that can induce fission for relativistic heavy-ion col- 
lisions. The predominant mechanism is nuclear fission. This is the case when the per- 
pendicular distance between ion centers as they pass each other, the impact parameter 6, 
is less than the sum of the ions’ radii. Another mechanism is due to the electromagnetic 
interactions between the ions. This is the case when b becomes greater than the sum 
of the radii of the two ions and it is called Coulomb fission (ref. 7). The objective of 
this study is to calculate the Coulomb fission cross sections for 208 Pb and 209 Bi ions at 



158 A GeV interacting with a stationary 208 Pb target using the full Weizsacker- Williams 
(WW) method and an approximate WW method. We will verify the Coulomb fission and 
the electromagnetic absorption cross section calculations of Abreu et al. (ref. 10) based 
on the approximate WW method that they used. We will be using their notations and 
conventions. 


2. The Equivalent Photon Approximation 


The Equivalent Photon Approximation, also known as the Weizsacker- Williams Method 
of Virtual Quanta (ref. 11), is a classical computational scheme based upon a plane wave 
approximation to radiation that is Lorentz contracted in the direction of motion and con- 
centrated normal to that direction. The original idea came from Fermi in 1924 (ref. 12 ). 
It was then extended independently by Williams in 1934 (ref. 13) and Weizsacker in 1934 
(ref. 14). Consider the target ion, 5, from the reference frame of the projectile ion, A. 
The stationary target has an apparent velocity toward the projectile ion. As the speed 
of closure approaches that of light, the target’s now electromagnetic (motional induced 
magnetic) field can be modeled as plane wave radiation or as an equivalent swarm of 
virtual photons. We view the collision in cylindrical coordinates centered on the target 
where the projectile trajectory is parallel to the axis but offset by the impact parameter, 
b. The target ion’s electric field is Lorentz contracted in the longitudinal direction along 
the axis and concentrated with respect to the circular radial direction. Accordingly, the 
photon radiation from the target can be replaced by two radiation pulses, Pi and P 2 . 
Pulse Pi has the concentrated electric field perpendicular to the motion and pulse P 2 has 
the longitudinal electric field. 

The Weizsacker- Williams method uses a function, N(E ), that represents the number 
of virtual photons with energy E in the radiation pulses per unit energy. This number 
spectrum can be written as an integral over the impact parameter plane as 


N(E) = 



N(E , b)2irbdb 


(i) 


where b m i n is the minimum impact parameter. The argument in equation (1) can be 
written in terms of the contributions from the two pulses Pi and P 2 for a point ion (ref. 
11 ) as 


N(E,b) 
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( 2 ) 


where a = ^ (the QED fine structure constant), £ = (3 = 7 is the Lorentz 

contraction factor, Z is the number of positive charges in the ion, and K n (£) are modified 
Bessel functions of the second kind of order n. Equation (2) gives the number of virtual 
photons (quanta) with energy E = hw at transverse position b from the target ion per 
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unit energy per unit area as seen by the projectile ion traveling at speed v. Note that 
as the speed v — » c, then ^ 0. In this limit the second term in equation (2) can be 

dropped. 

The integration of equation (1) can be performed by using the modified Bessel differ- 
ential equation with relations for the first derivatives of the modified Bessel functions of 
the second kind. The resulting expression becomes 

2rvZ 2 f r 2 /? 2 1 

N(E) = Jpz jrfoMATiM - ~Y~ [KiM 2 - A-„M 2 ] | (3) 


where x = $ b ™ in = . Notice in equation (3) that the contribution from the second 

pulse P 2 is multiplied by a factor of 4 *, which is hidden in the (3 2 = 1 — ^ coefficient. 
When (3 — > 1 as ^ c, the contribution from the second pulse vanishes. 

The process we are interested in is the Coulomb fission cross section of the projectile 
ion A, which breaks up as it passes the target ion B in the peripheral collision. This total 
cross section can be written using the notation of reference 10 as 




N B (E)<rfl (E)dE 


(4) 


where a A ( E ) represents the microscopic photofission cross section, C indicates Coulomb, 
/ fission, and 7 a single photon process. Equivalently, equation (1) can be substituted 
into equation (4) to find 


POO pEmax 

a A c f - / 2Mb / N b (E, b)a 7 A f (E)dE (5) 

” bmin ** Emin 

which gives another form for the total Coulomb fission cross section. 

3. The Photofission Cross Section 

The photofission cross sections cr 1 / (E) for A = 208 lead and A = 209 bismuth pro- 
jectile ions were constructed according to reference 10. Graphical functions for these 
photofission cross sections in the photon energy range from 20 MeV to 240 MeV were 
found in reference 15. However, the overall semilog graphical functions for the energy 
range from 100 MeV to 1000 MeV were found in reference 16. For calculations, we took 
photofission points every 20 MeV in photon energy E. These functions are represented 
by the semilog point graphs in figure 1 for 208 Pb and in figure 2 for 209 Bi. 

We extrapolated the photofission graphs by linear extensions of the semilog points to 
2000 MeV according to reference 10. The extrapolations are designed to capture the total 
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cross section results found in reference 10 when the approximate WW method of equation 
(9) is used. By reproducing their results we established photofission cross sections similar 
to those found in reference 10. To see how sensitive the total cross section calculation 
is to the particular extrapolation used we considered three cases with point indicators in 
comparison with reference 10: circle (high), square (similar), and diamond (low) as shown 
in figures 1 and 2. 

Figure 1 shows our points for the photofission cross section of 208 Pb with three extrap- 
olations beyond 1000 MeV. Figure 2 shows our points for the photofission cross section of 
209 Bi with three extrapolations beyond 1000 MeV. The middle extrapolations with square 
or similar point indicators result in Coulomb fission cross sections that compare well with 
those reported in reference 10. 

4. Calculation of the Coulomb Fission Cross Section 


The calculations of the Coulomb fission cross sections were done using a simple ex- 
tended trapezoidal integration scheme. The energy variable E was integrated from 20 
MeV to 2000 MeV. The number of points used for the microscopic photofission cross 
section was 100. The energy interval was 20 MeV to accommodate this range. The mimi- 
mum impact parameter for both collisions, Pb-Pb and Bi-Pb, was chosen to be b m i n =15 
fm to correspond with the value used by reference 10. Note the maximum photon energy 
occurs at b min and is given by E max = ^£. 

We calculated cross sections using the full WW method by substituting equation 
(3) into equation (4) and integrating numerically as described with a standard modified 
Bessel function routine. We calculated cross sections using the approximate WW method 
by dropping the second term in equation (2). This is the contribution from the radiation 
pulse P 2 that is associated with the longitudinal electric field. This is equivalent to 
assuming v = c or that (3 = 1 . The resulting expression contains the modified Bessel 
function Afi(£) 2 . According to the asymptotic behavior of this function, the approximate 
expression can be split between a low photon energy and high photon energy expression. 
These energy approximations (ref. 10) are represented by 


N B (E,b)~ 


Z B 2 a 

7T 2 (3 2 b 2 E 


(6) 


for the low energy approximation E <C and 


N B (E,b)c± 


Z B 2 a 

2^b eXP 



(7) 


for the high energy approximation E 
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Following reference 10 we used the low energy approximation of equation (6) and set 
Nb(E , fc) = 0 for energies E > That is, we only integrate up to the cutoff photon 

energy Substituting equation (6) into equation (5), this approximation to the full WW 
method gives the Coulomb fission cross section as 


c f 

° A = 
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Intergrating over the impact parameter b results in 


a C J = f EmaX aZf(E) In ( lEE) 
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( 8 ) 

( 9 ) 


Equation (9) was used to do the approximate WW numerical calculations corresponding 
to those of reference 10. 

Table 1 shows a comparison between our results and those of reference 10 for 208 Pb and 
209 Bi beam ions interacting with a 208 Pb target. Column 1 gives the projectile ion; column 
2 gives the results of reference 10. Our approximate calculations using equation (9) are 
shown in column 5. For 208 Pb the calculation based on the similar type extrapolation of 
figure 1 gives 380 mb, which is comparable to the Coulomb fission cross section calculated 
by reference 10. For 209 Bi the calculation based on the similar type extrapolation of figure 
2 gives 450 mb, which is comparable to the value calculated in reference 10. Calculations 
using equation (9) based on the extrapolations represented by circle and diamond are, 
respectively, higher and lower, when compared with those of reference 10. Convergence 
was checked by varying the number of integration points. 

The experimental Coulomb fission cross section of 649 mb, shown in table 1, is our 
calculation based on reference 10. They calculated an expected yield of Coulomb fission 
events per incident 208 Pb ion of 0.9 x 10~ 2 . The reference 10 calculation included the 
contributions of the isotopes 208 Pb, 207 Pb, and ^Pb, which were integrated through 
the 12-mm thickness of the target. It was assumed that all the isotopes had the same 
Coulomb cross section of 380 mb. From the expected yield it is necessary to subtract off 
an 18-percent estimated correction for nuclear reinteraction of the fission fragments inside 
the target. However, the observed yield of Coulomb fission events per incident 208 Pb ion 
was (1.26 d= 0.16) x 10 2 (from NA50 experiment at CERN SPS). Adjusting the observed 
yield number to before the 18-percent correction for nuclear reinteraction would be about 
1.536 x 10 2 , which is associated with an experimental cross section. This experimental 
Coulomb fission cross section becomes (^p)380 mb « 649 mb. 
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Tab 

le 1 . Comparison of Coulom' 

3 Fission Cross Sections 

Ion 

Ref. 10 

Based 

Figs. 1 & 2 

Approx. WW, mb 

Full WW, mb 

A 

Calc., mb 

on Exp., mb 

^J(E) 

Eq. (9) 

Eq. (3) in Eq. (4) 

208 p b 



Circle 

395 

331 

1 mug 

380 

649 

Square 

380 

315 




Diamond 

361 

296 

IBUfii 



Circle 

454 

372 

209 Bi 

450 


Square 

450 

369 

IHgg 



Diamond 

448 

368 


According to reference 10 the expected yield of Coulomb fission events per incident 
2 °8pb ion after subtraction of an 18-percent correction for fission fragments is about 
0.75 x 10 2 , which is about 40 percent lower than the observed yield. Based on our 
more accurate, full WW calculations for the Coulomb fission cross section for the 208 Pb 
projectile of 315 mb shown in table 1, the revised calculated yield becomes (|^j)(0.9 x 
10~” 2 )(1.00 — 0.18) « 0.61 x 10~ 2 . This is now about 52 percent below the observed yield 
of 1.26 x 10- 2 . 

5. Calculation of the Electromagnetic Absorption Cross Section 

The electromagnetic absorption (abs) cross section is calculated by replacing the 
photofission cross section a^(E) with the photoabsorption cross section a^ bs (E) in the 
previous equations. We constructed a rough log-log point graph of this absorption cross 
section for ^Pb from a graph found in reference 17 as shown in figure 3. Table 2 shows 
our results in comparison with reference 10. Column 3 data compare well with that of 
column 2, which are also based on equation (9). The full WW method shown in column 
4 results in a lower cross section. 


Table 2. Electromagnetic Absorption Cross Section 


Ion 

A 

Ref. 10 
Calc., b 

Approx. WW, b 
Eq. (9) 

Full WW, b 
Eq. (3) in Eq. (4) 

2 08pb 

50 

49 

44 


Integration of equations (4) and (9) was done over five photon energy ranges due to 
the form of our log-log graph for the absorption cross section: from 6 MeV to 20 MeV 
with an energy increment of 1 MeV; from 20 MeV to 200 MeV with an energy increment 
of 10 MeV; from 200 MeV to 300 MeV with an energy increment of 20 MeV; from 300 
MeV to 1000 MeV with an energy increment of 100 MeV; and from 1000 MeV to 2000 
MeV with an energy increment of 200 MeV. We used the same trapezoidal integration 
scheme and convergence check as that discussed in section 4. 
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6. Discussion 


The reason that the full WW method results in smaller Coulomb fission cross section 
compared with the approximate WW method is due to the slight difference in the respec- 
tive WW number spectrums as shown in figure 4. The N(E) s match at low values and 
high values of E , but through the central energy range the full WW number spectrum is 
offset lower than the approximate WW number spectrum. 

The calculated expected yield of Coulomb fission events per incident 208 Pb projectile 
ion was reported to be 40 percent lower than the experimentally observed yield in refer- 
ence 10. This was based upon their calculated Coulomb fission cross section of 380 mb. 
However, with our more accurate cross section of 315 mb, the discrepancy becomes worse: 
52 percent lower for the calculated yield versus the observed yield. For an experimentally 
based effective cross section of about 649 mb, this implies that the physics of the Pb-Pb 
interaction is still not understood. 

The key to application of the WW method is the theoretical or experimentally based 
photonuclear cross section. Our photofission cross sections of figures 1 and 2 are re- 
constructions from inferred cross sections based on experimental electron-induced fission 
cross sections (ref. 16). After 1000 MeV we used linear semilog extrapolations up to 2000 
MeV. Why the 208 Pb graph is rising and the 209 Bi graph is decreasing at 1000 MeV would 
reflect the nature of the fission for these nuclei. 

The choice of b m i n has an important effect on the resulting cross section. Smaller num- 
bers for b m i n result in larger cross sections. In order to avoid inducing strong interactions 
bmin = Ra + Rb- We used 15 fm for the minimum impact parameter, the value used 
in reference 10. We also looked at values of b m i n based on other methods (refs. 4, 18, 
and 19). These methods resulted in about the same discrepancy as reported previously. 
Calculations using b m i n according to the Wood-Saxon method (refs. 20 and 21) resulted 
in larger cross sections, however, not large enough to remove the discrepancy. 

The interference between nuclear and Coulomb forces is expected to be small because 
of the different distance behavior of the forces: the Coulomb force being weak and long 
range and the nuclear force being strong and short range. (This is discussed in ref. 18, 
where they show that the interference is small.) Thus, the question over the discrepancy 
between the calculated expected yield per Coulomb fission event and the experimentally 
observed yield remains open. However, because of previous excellent agreement with 
theory and experiment, we suspect that the solution of the disagreement lies with the 
complicated and approximate procedure used to extract the experimental cross sections. 


7 



References 


1. Fasso, A.; Ferrari, A.; Ranft, J.; and Sala, P.: FLUKA-99. Official web site of 
FLUKA code: http://fluka.web.cem.ch/fluka/ Accessed October, 4, 2001. 

2. Ferrari, A.; Sala, P.R.: The Physics of High Energy Reactions. Workshop on Nuclear 
Reaction Data and Nuclear Reactors. International Center for Theoretical Physics 
( Miramare- Trieste, Italy, 1996). 

3. Waters, L.S.: MCNPX. Official web site of MCNPX code: http://mcnpx.lanl.gov/ 
Accessed April 20, 2000. 

4. Bertulani, C.; and Baur, G.: Electromagnetic Process in Relativistic Heavy Ion 
Collisions. Phys. Rep., vol. 163, 1988, pp. 229-408. 

5. Norbury, J.W.; and Baur, G.: Explanation of Recent Observations of Very Large 
Electromagnetic Dissociation Cross Sections II. Higher Order Corrections. Phys. 
Rev. C, vol. 48, 1993, pp. 1915-1918. 

6. Wheeler, R.; and Norbury, J.W.: Higher Order Corrections to Coulomb Fission. 
Phy. Rev. C ' vol. 51, 1995, pp. 1566-1567. 

7. Norbury, J.W.: Relativistic Coulomb Fission, Phys. Rev. C, vol. 43, 1991, pp. 
368-370. 

8. Norbury, J.W.; Cucinotta, F.A.; Townsend, L.W.; and Badavi, F.F.: Parameterized 
Cross Sections for Coulomb Dissociation in Heavy-Ion Collisions. Nucl. Instr. Meth. 
Phys. Res. B , vol. 31, 1988, pp. 535-537. 

9. Wilson, J.W.; Tripathi, R.K.; Cucinotta, F.A.; Shinn, J.L.; Badavi, F.F.; Chun, 
S.Y.; Norbury, J.W.; Zeitlin, C.J.; Heilbronn, L.; and Miller, J.: NUCFRG2: An 
Evaluation of The Semiempirical Nuclear Fragmentation Database. NASA TP-3533, 
1995. 

10. Abreu, M. C.; et al.: Observation of Fission in Pb-Pb Interactions at 158 A GeV, 
Phys. Rev. C , vol. 59, 1999, pp. 876-883. 

11. Jackson, J. D.: Classical Electrodynamics. Third ed., John Wiley & Sons, New 
York, 1999. 

12. Fermi, E.: Uber die Theorie des Stofles Zwischen Atomen und Elektrisch Geladen 
Teilchen. Z. Physik , vol. 29, 1924, pp. 315-327. 

13. Williams, E.J.: Nature of the High Energy Particles of Penetrating Radiation and 
Status of Ionization and Radiation Formula, Phy. Rev. Let., vol. 45, no. 11, 1934, 
pp. 729-730. 



14. Weizsacker, C.F.: Ausstrahlung bei Stofien Schneller Elektronen. Z . Physik , vol. 
88, 1934, pp. 612-625. 

15. Arruda-Neto, J.D.T.; Sugawara, M.; Miyase, H.; Kobayashi, T.; Tamae, T.; Abe, 
K.; Nomura, M.; Matsuyama, H.; Kawahara, H.; and Namai, K.: Electrofission 
208 Pb in the Intermediate Energy Region. Phys. Rev. C. vol. 41, no. 1, 1990, pp. 
354-357. 

16. Moretto, L. G.; Gatti, R.C.; Thompson, S.G.; Routti, J.T.; Heisenberg, J.H.; Mid- 
dleman, L.M.; Yearian, M.R.; and Hofstadter, R.: Electron- and Bremsstrahlung- 
Induced Fission of Heavy and Medium-Heavy Nuclei. Phys. Rev., vol. 179, 1969, 
pp. 1176-1187. 

17. Krauss, F.; Greiner, M.; and Soff, G.: Photon and Gluon Induced Processes in 
Relativistic Heavy-Ion Collisions. Prog. Part. Nucl. Phys., vol. 39, 1997, pp. 
503-564. 

18. Benesh, C.; Cook, B.; and Vary, J.: Single Nucleon Removal in Relativistic Nuclear 
Collisions. Phys. Rev. C ', vol. 40, 1989, pp. 1198-1206. 

19. Kox S.; Gamp, A.; Perrin, C.; Arvieux, J.; Bertholet, R.; Bruandet, J.F.; Buenerd, 
M.; Cherkaoui, R.; Cole, A.J.; El-Masri, Y.; Longequeue, N.; Menet, J.; Merchez, F.; 
and Viano, J.B.: Trends of Total Reaction Cross Sections for Heavy Ion Collisions 
in the Intermediate Energy Range. Phys. Rev. C , vol. 35, 1987, pp. 1678-1691. 

20. John, S.; Townsend, L.; Wilson, J.; and Tripathi, R.: Geometric Model From Mi- 
croscopic Theory for Nuclear Absorption. NASA TP-3324, 1993. 

21. Krane, K.: Introductory Nuclear Physics. John Wiley & Sons, New York, 1988. 


9 



10 1000 1500 2000 

Photon Energy, MeV 

*ure 1. Photofission cross section for 208 Pb. 
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Figure 3. Photoabsorption cross section for 208 Pb. 
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Figure 4. Comparison of the full (solid line) and approximate (dashed line) WW number spectrums per unit 
energy. 
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